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I. INTRODUCTION 



The SELEX collaboration at Fermilab recently reported evidence of five resonances that 
may possibly be identified with doubly charmed baryon states Tentatively the states have 
been interpreted as ccd+(3443), ccci+(3520), ccw ++ (3460), ccw++(3541) and ccw ++ (3780). 
Subsequently the ccd + (3520) state has been confirmed in two different decay modes (H+ — > 
A~^K~7T + ; — > pD + K~) at a mass of 3518. 7±f .7 MeV with an average lifetime less than 33 
fs. Although these findings need to be confirmed by other experiments and larger statistical 
samples, they have triggered a renewed theoretical interest in doubly heavy baryon systems. 

Doubly heavy baryons have been studied with several methods, mostly non-relativistic 
ootential models (for some reviews see P, Oh, but also relativistic models P| 
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been obtained from lattice calculations 




sum rules 

and in a chiral Lagrangian framework I 7|. M asses of the lowest lying resonances have 
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121 ] . Doubly heavy baryons are also 



suited to be studied in a QCD Effective Field Theory (EFT) framework. Indeed, they 
are characterized by at least two widely separated scales: the large heavy-quark masses, 
m, and the low momentum transfer between the heavy and the light quarks, which is of 
order Aqcd- If one assumes that the typical momentum transfer between the two heavy 
qual ,s is >a rg e r than A QCD , then a QQq ha ry0n is ve ry shnhar to a hound state of a heavy 
antiquark and a light quark. This has first been noted in [13J, where at leading order in 
AqcdA^ the hyperfine splitting of the doubly heavy baryon ground state has been related 
to the ground-state hyperfine splitting of the heavy-light meson. In 14| non-leptonic and 
semileptonic decays of doubly heavy baryons have been examined in the context of SU(3) 
flavour symmetry. After this original work no further step has been made in the direction of 
providing a systematic description of doubly heavy baryons in an EFT framework that fully 
combines the dynamics of the two heavy quarks with that one of the light one. Following 
some suggestions in Q|> with this work we attempt to make such further step. In particular, 
we identify the degrees of freedom and write the low-energy EFT Lagrangian that describes 
doubly heavy baryon systems in the heavy-quark sector, once the heavy-quark momentum- 



transfer scale has been integrated out. The framework is similar to that one developed in 

n 

the last years for heavy- quarkonium systems (for a review see |lfij). 

Baryons made of three heavy quarks QQQ have not been observed yet. Their relevance 
has been emphasized since long ago [17]. They would reveal a pure baryonic spectrum 
without light-quark complications and provide valuable insight into the quark confinement 
mechanism. Indeed, the three-quark static Wilson loop is intensively studied on the lattice 
Jisl lfll ] as a source of information about the baryon heavy-quark potential and the type of 
confining configurations In this work we will identify the degrees of freedom 

and write the low-energy EFT Lagrangian that describes heavy baryons made of three heavy 
quarks, once the heavy-quark momentum-transfer scale has been integrated out. We will 
express the leading-order and spin- dependent potentials in terms of Wilson loop amplitudes 
along the lines developed for heavy quarkonia in 23] ]. 

A recent review that also discusses the present status of the art, experimental and theo- 
retical, including lattice, for heavy baryons made with two or three heavy quarks is Ref. 2^ ]. 
We refer to it for a more complete bibliography on the subject. This work is partially based 
on (25]. We refer to it for details in some of the derivations. 

The paper is distributed as follows. In Sec. HTIwe introduce NRQCD for heavy baryons. In 
Sec. II I II we write the low-energy EFT for QQq baryons and calculate the hyperfine splitting 
of the ground state. In Sec. IIVI we write the low-energy EFT for QQQ baryons and give 
some exact non-perturbative expressions for the leading-order and spin-dependent potentials. 
Sec. El is devoted to the conclusions. Some technical details may be found in the appendices. 



II. NRQCD 

Non-relativistic QCD (NRQCD) is the EFT suitable to describe systems made of two or 
more heavy quarks. It is obtained from QCD by integrating out modes of energy of the order 
of the heavy-quark masses j^ . 

We are interested here only in the heavy-quark sector of the NRQCD Lagrangian. The 
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2-heavy-quark sector coincides with the Lagrangian of the Heavy Quark Effective Theory 
(HQET). Up to order 1/m 2 it reads: 

N Q r n2 



£ NR Q CD = £ Q t 



D 2 {h) <r-gB {h) [D;gE] , , (h) a ■ [D x , gE\ 
tD + - — + cV — + & D > + i& s — 2 



L 2m h ' -* 2m h ' ~ u 8m 2 h ' ~~ a 8m 2 h J Qh ' ^ 

where Nq is the number of heavy-quark flavours, Qh the Pauli spinor field that annihilates the 
quark of flavour h and mass m^, iD = id$ — gA°, iD = zV + gA, [D-, E] = D ■ E — E ■ D, 



for instance, in 



» » n „A >' 



[Dx,E] — DxE — ExD, E % — F i0 , B l = -e ijk F^ k /2 (e 123 = 1) and a = (tr u <r 2 , a 3 ) 



are the Pauli matrices. The matching coefficients c F , c D and c s may be found at one loop, 



271 ] . The Lagrangian (0) (with 0(l/m 3 ) terms included, but all matching 



coefficients set equa 



to 1) has been used to perform lattice calculations of the spectra of 

heavy baryons in 

At order 1 / m 2 the NRQCD Lagrangian relevant to describe baryons made of two or more 
heavy quarks exhibits also a 4-heavy-quark sector: 
Nq / jss jsv 

(In, r> . j. i tlf 



/'NRQCD 

QQ 



/ d ss d sv 
h ,fj? =1 \ m hm h > m h m h , 



d vs " d vv ^ \ 

+ ^Q JlL Q[T a Q h Ql,T a Q h/ + V Q{T a crQ h • Q^TVQ,, . (2) 

m h m h , m h m h > ^ ) 

a=l a=l / 

The matching coefficients dQ h Q , start getting contributions at order a 2 . They have been 
calculated to this order in appendix 1X1 

Six-quark operators contribute to heavy baryons made of three heavy quarks. They show 
up at order 1/m 5 . The corresponding matching coefficients start getting contributions at 
order a*- Hence, these operators are highly suppressed and will be neglected in the rest of 
the paper. 

In the following we will make a step further and construct the EFT suitable to describe 
baryons made of two (Sec. IIII|) and three (Sec. IIV|) heavy quarks once gluons of energy or 
momentum of the order of the momentum transfer between the heavy quarks have been 
integrated out. The procedure and the resulting EFT will be quite similar to that one 
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developed for heavy quarkonium in 



28, 



2911 . For this reason we will call the EFT with the 



same name: potential NRQCD (pNRQCD). 



III. PNRQCD FOR QQq BARYONS 

In this section we deal with baryons made of two heavy quarks Q\, Q2 (bb, be or cc) with 
masses mi and m 2 respectively and one massless quark q. The dynamics of these systems 
is expected to mix aspects typical of heavy quarkonium with aspects typical of heavy-light 
mesons. On the one hand, the interaction of the two heavy quarks is that one of a non- 
relativistic quark pair close to threshold moving with relative velocity v. It is, therefore, 
characterized by the energy scales: m ^> mv ^> mv 2 , where mv is the scale of the typical 
momentum transfer between the two heavy quarks (or of the inverse of their typical distance) 
and mv 2 is the typical binding energy. On the other hand, the energy scale that governs the 
interaction between the heavy quarks and the light one is Aqcd- Two different situations 
are possible. 

(A) If mv ^> Aqcd, at a scale /1 such that mv ^> fi ^> Aqcd the heavy-quark distance 
cannot be resolved. The Q1Q2 pair behaves like a point-like particle (sometimes also called 
diquark j^O]) in an antitriplet or sextet colour configuration. In the antitriplet configuration 
the two heavy quarks attract each other. The interaction of the antitriplet field with the 
ight quark is similar to that one of a heavy antiquark with a light quark in a D or B meson 
However, the spectrum is expected to be richer and more complex due to the internal 
excitations of the heavy-quark system. These include radial and spin excitations, but also 
colour excitations to sextet configurations. Considering that the scale /1 is perturbative, 
in this situation pNRQCD has the following degrees of freedom: light quarks, gluons of 
energy and momentum lower than mv (also called ultrasoft) and heavy quarks. To ensure 
that the gluons are of energy and momentum lower than mv , gluons appearing in vertices 
involving heavy-quark fields are multipole expanded in the relative distance r ~ l/(mv) 
between the two heavy quarks. This corresponds to expanding in rA QC D ~ ^QCD/( m v) 
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or r (mv 2 ) ~ v. In the situation mv ^> Aqcd, in principle one may further distinguish 
between the subcases mv 2 ^> Aqcd, mv 2 ~ Aqcd and Aqcd 3> mv 2 . In the first case, one 
may expand in Aqcd/(™^ 2 ) and disentangle the heavy-heavy dynamics, which is completely 
accessible to perturbation theory, from the heavy-light one. In general, excitations of the 
heavy-heavy system will dominate over excitations of the heavy-light system. In the latter 
case, the potential governing the heavy-heavy system gets non-perturbative contributions. 1 
Since the kinetic energy of the heavy quarks is smaller than Aqcd, the heavy-light dynamics 
dominates in this situation over the heavy-heavy one. At leading order in the mv 2 / 'Aqcd 
expansion the flavour symmetry typical of the HQET is restored. 

In the rest of this section, we will deal with the general situation mv ^> Aqcd, without 
assuming any special hierarchy between the scales mv 2 and Aq C d- To be definite, one may 
think that we work in the situation mv 2 ~ Aqcd- 

(B) If mv ~ Aqcd the distances between the three quarks are of the same magnitude. 
Hence, we cannot disentangle the heavy-quark pair dynamics from the light-quark one. 
Moreover, the potential between the two heavy quarks is non-perturbative. At the level 
of NRQCD, the system may be studied with lattice calculations. In this situation it seems 
unlikely that a simple diquark-light-quark picture holds. In general, from an EFT point 
of view it does not seem consistent to have a diquark-light-quark picture for the heavy- 
quarks-light-quark interaction, which implicitly assumes mv ^> Aqcd, and at the same time 
a confining potential binding the two heavy quarks, which requires mv ~ Aq C d, as so often 
done in potential models. 

In the following, we will work out pNRQCD in the situation (A). This situation is ex- 
pected to be appropriate for the description of at least doubly heavy baryons in the ground 
state. 

1 An analogous situation for the quarkonium system has been treated in |29| . 
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A. Lagrangian 



In this section we write the pNRQCD Lagrangian that describes heavy baryons of the 
type Q1Q2Q in the situation where the typical momentum transfer between the two heavy 
quarks is much larger than Aqcd- This corresponds to the case labeled (A) above. 

The number of allowed operators is reduced if we choose to have a manifestly gauge 
invariant Lagrangian. This may be obtained by projecting the Lagrangian on the heavy- 
heavy sector of the Fock space, by splitting the heavy-heavy fields into an antitriplet and 
sextet component (r = X\ — x 2 , R = {m\X\ + m 2 a^ 2 ) / (mi + 7712)), 



Qu{x u t)Q 2j {x 2 , t)~Y, T '( r ' *) ft + E S ^ r ' fl ' *) ^' 3 = 2 ' 3 ' ( 3 ) 



and by building the Lagrangian from these operators. The tensors T^- and S?- are defined 
in appendix IB II The Lagrangian is constrained to satisfy all the symmetries of QCD. 
In particular, it is symmetric under the exchange of the heavy quarks. Such symmetry 
transformation changes mi <-> m 2 and r to — r. The gluon fields are even, because multipole 
expanded around the centre-of-mass of the heavy-heavy system. For what concerns the 
heavy-quark fields, from Eq. (jSJ) it follows that T e is even, because ■ is odd under the 
exchange i <-> j, and S CT is odd, because S^- is even under the exchange i <-> j. 

The resulting Lagrangian £ p nrqcd = £ p nrqcd(-R, a ^ 0{l/m) in the 1/m expansion 
and at 0(r) in the multipole expansion is 
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1=1 



(7=1 



£pNRQCD — £ gluon 




pNRQCD' 



(4) 



with 




(5) 



a=l 
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(6) 
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,.(0,0) 
^pNRQCD 



d 3 rT f 



iD - \4 0) 



T + E f 



iD - K 



(0) 



(7) 



(o,i) 

pNRQCD 



— / 



TV .EE 



8 3 6 

EEE 

a=l ^=1 <t=1 



\ijk=l 



Vij'fc=l 



^1 ~ m 2 T/ (0,l) 

2m^ 



^ E TV • T - E sV • s > ( 8 ) 



a=l 



,(1,0) 
-pNRQCD 



(1,0) 



d 3 rT f 
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T + E f 
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_2m R 


2m r 


_2m R 


2m r 





TCTBY, 



EEE 



a=l £=1 <x=l 



E ) T ^ 



2mi 



+ 



2m 9 



gB a YT 



V, 



(1.0) 



1(7 BT jif 

o=l 

(1,0) 8 
£<X BE 



E 251531 E' 

ijk=l / 

+ -£ 

2mi 2m 2 



c>'cr 



c>'cr 



a=l 



2mi 
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cytr- 



2mi 

gB a T a T 

gB a T£ E 



+ 



c F 'cr 



2m 2 



+« E ^EEE 



a=l €=1 cr=l 



ajfc=l 



T/(1,0) 

1 TL-BT 
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T/(1,0) 
1 ^ELBE 
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- E^^Tl E^-^T^ 

\ijfc=l / 
9 \ 8 

E TtL - -gB a T§T 



a=l 



2 \ 



(9) 



a=l 



where m R = mi + m 2l m r = mim 2 / (mi + m 2 ), (X^ 1 ) is the Pauli matrix acting on the heavy 
quark h, iD R = iV R + gA, L r = r x (-iV r ), T = {T l ,T 2 ,T 3 ), E = (E 1 , E 2 , E 6 ), the 
gauge fields in the covariant derivatives acting on the antitriplet and sextet are understood 
in the antitriplet and sextet representation respectively, T| and T 6 a have been defined in 



appendix IB 31 and all gluon fields are evaluated in (R, t) . The coefficients Cp and Cp are 
the Wilson coefficients of NRQCD introduced in Sec. |H] The functions V are the Wilson 
coefficients of pNRQCD for doubly heavy baryons. They encode the contributions coming 
from gluons of energy or momentum of order mv, which have been integrated out. They are 
non-analytic functions of r. As we will discuss in the next section, at tree level they are 



V, 



(o,i) 



Tr-EE 



(0,1) 



(1.0) 



TO BY. 



V 



TrET 



(1.0) 
T(T BT 



v 



(0,1) 



v 



Et--.EE 



(1,0) 
E<T BE 
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Vn 



(1,0) 



TLBS 



V 



(1,0) 



TL-BT 



V 



(1,0) 



ELBE 



while Vj^ and get the first non- vanishing contribution at order a s . In Eqs. (JHJ) and (JSJ) we 
have displayed only the operators that have a non- vanishing tree-level matching coefficient. 
The coefficients in front of the D 2 R and operators in (jHJ) are equal to 1, due to Poincare 



3l|. We observe that 



invariance or dynamical considerations similar to those developed in 
in the case m\ ^ m 2 , electric dipole transitions between antitriplet states induced by the 
term 

mi — m 2 

a=l 



are allowed j], 3^ . 



2m 



R 



^T ] r ■ gE a T?T 



The power counting of the Lagrangian (j3J) in the centre-of-mass frame goes as follows: 



mv, r ~ l/(mv), D R ~ Aqcd, mv 2 , V?% ~ 



mv 2 and E, B 



A^cd, (mv 2 ) 2 . The 



power counting is not unique, because the scales mv 2 and Aqcd are still entangled in the 
dynamics. The Lagrangian at leading order reads 



^-pNRQCD 



rf 3 rT f 
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iD 



2m r 



V 



(0) 



iD 



2m r 



(o) 



(12) 



a=l 



/=1 
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B. Matching 



The matching from NRQCD to pNRQCD is, in general, performed by calculating Green 
functions in the two theories and imposing that they are equal order by order in the inverse 
of the mass and in the multipole expansion. Since we are working in the situation where 
the typical momentum transfer between the heavy quarks is larger than Aqcd ; we can, in 
addition, perform the matching order by order in a s . 

If we aim at calculating the matching at tree level a convenient approach consists in 
projecting the NRQCD Hamiltonian on the two-quark Fock space spanned by 

r 3 

/ d 3 Xl d 3 x 2 ^^ iQ2 (x 1 ,x 2 )Q\\x 1 )Qi\x 2 )\0), (13) 
J ij=i 

where |0) is the Fock subspace containing no heavy quarks but an arbitrary number of 
ultrasoft gluons and light quarks and &q 1 q 2 (x 1 , x 2 ) is a 3 (g> 3 tensor in colour space and 
a 2 <g> 2 tensor in spin space. This is similar to what is done in ^j]. After projection, all 
gluon fields are multipole expanded in r. In order to make gauge invariance explicit at 
the Lagrangian level, it is useful to decompose ^q 1 q 2 (x\,x 2 ) into a field T(r, R,t), which 
transforms like a colour antitriplet, and a field H(r, R,t), which transforms like a colour 
sextet: 

3 

^QiQa^l' 3 ^) = E <l>ii'(xi, R,t)(j) jf (x 2 ,R,t) 

'3 6 \ 



,e=i (t=i 



where 



(y, x,t) = P exp fig J ds {y - x) ■ A(x + (y - x)s, t) \ . 



(15) 



P stands for path ordering. At leading order in the coupling constant, <f>ij(xi,R,t) = 5ij 
and 

3 6 

$g iQ2 ( Xl , x 2 , t) « e T '( r ' R > *) % + E SCT ( r ' K ' *) ^- ( 16 ) 



(7=1 
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After projecting on (|13|). one obtains the Lagrangian ©-© with the matching conditions 

(HOD- 

As an example, let us consider the calculation that leads to the term 

* W» - J *r\ g T» + • T (17) 

in the pNRQCD Lagrangian (see Eq. (JHJ)). We start from the NRQCD term 

^nrqcd = Q|cy — Qi + Q^cy — Q 2 • (18) 

Projecting onto we obtain in the antitriplet-antitriplet sector 

8 3 



«/ ^, — 1 DDI X 7 1 \ 1 



o=l £l'ijk=l 



+T^.^-.^^T^). (19) 



Using the definition (|Blj) . we have 
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E T^T£ = ~f = (20) 

ijk=l 

±.ij J -jk±-ik - 2 - 2 , 

and eventually end up with Eq. (JT7J). This fixes V^} BT = 1 at leading order. Note that 
Eq. ()17j) differs by a factor 1/2 from Eqs. (9) and (10) in [Lj], which seem to miss the correct 
colour normalization of the antitriplet states. 

One may ask what happens to Vj^f BT beyond tree level. Order a s corrections may only 
come from one-gluon corrections to the NRQCD vertex of Eq. (|18)l . because all other spin- 
dependent operators in NRQCD contribute to higher orders in 1/m. One-loop corrections 
to the external (transverse) gluon or to a quark line or involving a gluon attached to the 
external gluon and to the quark line coupled to it vanish in dimensional regularization, once 
we have expanded in the external energies. Gluons attached to a quark line are longitu- 
dinal. It is convenient to use the Coulomb gauge. In Coulomb gauge, longitudinal gluons 
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exchanged between different quark lines cancel in the matching with equal contributions 
from the pNRQCD side. Finally, longitudinal gluons attached to the external gluon line and 
a heavy quark line not coupled to it contribute to higher-order operators, since the three- 
gluon vertex is proportional to the external energies. We conclude that V^} BT does not 
get contributions at one loop, so that V^^ BT = 1 + 0(a^). Similar considerations hold for 
y (1,0) and T/ (1 ' 0) 

V SCT BE dIlu Tff BS 1 I - 

The perturbative matching of the static potentials and goes as follows (see 2£| 
for the quarkonium case). In NRQCD we compute static Green functions, whose initial and 
final states overlap with the antitriplet and sextet fields in pNRQCD. Since we work order 
by order in a s , it is not necessary for the Green functions to be gauge invariant. A possible 
choice is 

3 

IZ = E (0\Mt j Qi(R,xi,T/2)Q j (R,x 2 ,T/2)Ml f QUR,y 1 ,-T/2)Q] f (R,y 2 ,^ 

(22) 

(1) if M =T, M£ = T£, u,i; = 1,2,3, 

(2) if M = E, M^=^, u,v = l,2,....,6, 

where 

Q(R, x, t) = (j)(R, x, t)Q(x, t), (23) 

and 4>{R, x,t) has been defined in Eq. (jl5j) . Integrating out the static-quark fields from 1^ 
we obtain 

IZ = S 3 (x l - yi )5 3 (x 2 - y 2 )<0|(W^r|0> (24) 
with Wqq diagrammatically represented in Fig. [T] and explicitly given by 

3 

(w^ Q r = p J2 ^•0«'(fl,«i > r/2)0 w (r/2,-r/2,! Ba )0*'*(a5i,-R,-r/2) 

ijkni' j'k'n' =1 (25) 

x <f)jj>(R, x 2 , T/2)<f> jlnl (T/2, -T/2, x 2 )<f> nln (x 2 , R, -T/2) ML- 
In the large T limit, the Green functions 1™ and I™ are reduced to the antitriplet and 
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FIG. 1: Static Wilson loop with edges x 1 = (x 1 ,T/2), x 2 = (x 2 ,T/2), yi = (x 1 ,-T/2), y 2 = 
(x 2 ,-T/2) and insertions of the tensors My and Mj>j> in X = (R,T/2) and Y = (R,-T/2) 
respectively. 

sextet propagators of pNRQCD respectively. If we neglect subleading loop corrections to the 
pNRQCD side of the matching, we obtain: 



lim (0|(W&r|0) = lim Z M (r)e*p (-iV^^T) 

3 

x(0| ^jMT/2,-T/2,R)<l> jn (T/2,-T/2,R)ML\0), 

ijkn=l 



(26) 



wixvrv Z\t is a normalization factor. At order n s wo owl up with the well-known result j^: 

(27) 
(28) 



^(0 = -^, 
3 \r\ 

1 « s 



3 r 



The antitriplet channel is attractive, the sextet one repulsive. 



C. Hyperfine Splitting 

In the dynamical situation considered here (case (A) of Sec. IIII[) . a doubly heavy baryon 
is mainly a bound state of a heavy quark (or antiquark) pair in an antitriplet (or triplet) 
configuration and a light quark (or antiquark). Sextet field configurations show up in loops 
with ultrasoft gluons. Their contribution is suppressed either in the multipole expansion or 
in 1/m (see Eqs. (jHJ and 0). The leading-order pNRQCD Lagrangian is shown in Eq. (jf 2j) . 
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It does not depend on the spin of the heavy quarks. As a consequence, QQq baryons will 
appear in degenerate multiplets of the total spin Sqq^ = Sqq + Si, where S; is the spin of 
the light degrees of freedom and Sqq of the heavy-quark pair. This symmetry is similar 
to the spin symmetry of the HQET. Differently from the HQET, however, the pNRQCD 
Lagrangian depends at leading order on the heavy-quark flavour. This is a consequence of 
the fact that we cannot, in general, neglect the kinetic energy. 2 

We will consider in this section the S'-wave ground state of a doubly heavy baryon made 
of two identical heavy quarks Q of mass itlq. In this case, since an (anti)triplet state is 
antisymmetric in colour, due to the Fermi statistics, the two heavy quarks are allowed only 
in a spin 1 (symmetric) state. In the standard notation, the lowest energy states for QQu or 
QQd are called Sqq (Hqq) for spin 1/2 (3/2), and for QQs, Qqq (Qqq). Since the heavy- 
quark pair spin is fixed, the hyperfine splitting may only originate from spin-dependent 
couplings of the heavy quarks with the light one. The leading-order operator (in a Aqcd/?ti 
expansion) is given by Eq. (JT7j) . We will derive a simple formula that relates at leading order 
in the Aqcd/^ expansion the hyperfine splitting of a QQq doubly heavy baryon ground 
state with the hyperfine splitting of a Qq heavy-light meson ground state. The framework 
will be that one of pNRQCD, developed in the previous sections. The calculation will be 
similar to that one of Ref . . 

Let us consider, first, the case of a heavy-light meson Qq. The heavy antiquark may be 
described by a two-component field Q c = icr 2 Q*, where Q is the Pauli spinor that annihilates 
the heavy quark. We rename Q\ = Q + and Q\ = Q- since <2±|0) = \S% — ±1/2). Sq is 
the spin of the heavy antiquark and Sq q the total spin of the meson. The leading-order 
HQET Lagrangian does not contain spin-interaction terms, therefore, states that differ only 
in the spin quantum numbers are degenerate. In particular, this happens for the three 
lowest Sq q = 1 states (SI = 1, 0, —1), which we denote by \Pq), and for the lowest Sq q = 
state, which we denote by \Pq). An expression for these states that makes explicit their 

2 An exception may be the special case Aqcd 3> mv 2 . 
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heavy (anti)quark field content is given in appendix O The Hamiltonian responsible for the 
leading contribution to the hyperfine separation is 

^HQET = -Cf / d*R j2Qi(R) a ' 9 f {R)Ti Qc(R) 

J a=l Q 

8 

J d 3 R [{QllfQ+ ~ Q-TiQ-) 9B 3a + i(QlTiQ + - QlTiQ_) gB 2a 



2m Q 



+(QlT$Q- + QlTiQ + )gB la ], (29) 



where, after the last equality, we have dropped the explicit coordinate dependence of the 
fields. From Eq. ()29|) and Eqs. ()C1|) - (|C4|) it is straightforward to derive: 

1Q) r 8 
(P*\5H HQET \P^) - (P Q \SH mET \P Q ) = -2^- d 3 R(Sf = l/2\J2gB 3a n\Sf = l/2).(30) 

m Q J a=l 

In the case of a doubly heavy baryon QQq we proceed in a similar way. The triplet field 
T is a 2 (g) 2 tensor in spin space, which may be decomposed as 2 <g> 2 = 1 © 3, i.e. in a scalar 
component, T^ s \ and a vector one, T^: 

T y (r, R, t) = T^(r, R,t) + g (^) T^ k (r, R, t), i, j = 1, 2. (31) 

The indices ij refer to the spin space. Note that the matrices (ia 2 /V2)ij and (ia k a 2 / 'v2)ij 
are respectively antisymmetric and symmetric in ij. It is convenient to rewrite the fields 
T^ k as 

r = T^ 3 and T± = £^ + ^2! (32) 

v2 

since Tq|0) = |Sqq = 0) and 7±|0) = \Sqq = ±1). As we argued above, the leading-order 
pNRQCD Lagrangian describing doubly heavy baryons does not contain spin-interaction 
terms, therefore, states that differ only in the spin quantum numbers are degenerate. In 
particular, this happens for the four lowest Sqq q = 3/2 states (SqQ q = ±3/2, ±1/2), which 
we denote by |Sqq), and for the two lowest SQQ q = 1/2 states (Sqq ? = ±1/2), which we 
denote by \Eqq). An explicit expression of these states in terms of heavy (anti)triplet fields 
is given in appendix O The Hamiltonian responsible for the leading contribution to the 



15 



hyperfine separation is 



1Q) 



8K 



pNRQCD 



2m Q 

2rn Q 

J.Q) 
2m Q 



d 3 R / d 3 r K (1,0) 



T(T-BT\ 



gB a (R)T^ T(r, R) 



a=l 



d 3 R I d 3 r Vj^l^ r 



7<TB/' 1 'J E/ ^ ^ 

a=l ikj=l 



gB ka T a T (V)j 



d 3 R / d 3 r V, 



(1,0) 



TCT-BT 



a=l 



(TlT§T + - TlT$T_) gB 



] rpar 



>3 a 



+ ^(-TjT3 a T + 2jT|T + + tIT-£T - T$T§T_) gB 2a 
V2 

+ -^(TlT^T + T^TgT + + T^T^Tq + TIt$T_) gB la 
V2 

where the second equality follows from Eq. (j31|) and the third one from Eq. 
Eq. and Eqs. ()C5|) - ()C10|) it is straightforward to derive: 



(33) 
From 



( S QQ I ^pNRQCD | Sqq) - (Eqq] \5 H p nRQCd\Zqq} ~ 

r (Q) r 8 r 

~ 3 2^J d3R{S < = V2|g^ 3a 75lS? = 1/2) J d 3 r^ QQ (r)V^ 0) BT (r)^QQ(rUM) 

where lpqq is the ground-state eigenfunction of — Vj?/(2m r ) + V^°\ At NLO V^jj BT = 1, 
therefore f *r V ' QQ (r) V&VW VQQ (r) = l + This result crucially depends on the 

fact that we have multipole expanded the gluon fields. As a consequence, B does not depend 
on r and the magnetic dipole transition term (differently from the electric one oc r ■ gE) 
does not exhibit any explicit dependence on r. Comparing Eq. (|3U|) with Eq. (|34|) we obtain 
(Ms and Mp are the baryon and meson masses respectively) 



M=* 



Me 



QQ 



3m Q / cy 
4m Q c 



(Q) 

£_ ( M P * - Mp , 

(Q0 V Q' Q 



i + oK 2 , Aqcd Aqcd ^ 



(35) 



m Q niQi J 

Up to a factor 1/2, the formula is the one derived in In the previous section, the 

origin of the discrepancy has been traced back to a missing colour normalization factor 
in the spin antitriplet interaction term (JT7j) (and (|33|l). On the other hand, the relation 



M=* - M-, HS 

^QQ ^QQ 



3m Q > 



Mp. — Mp ) has been derived since long in non-relativistic po- 

Attiq \ Q' Q J 

tential models. Surprisingly the discrepancy between this formula and the formula in [13j 
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has to the best of our knowledge never been noticed before in the literature. 3 Even more 
surprisingly some of the literature has explicitly claimed agreement between the potential 
model prediction and the formula in 2|! 

From we read that M D * - M D = 142.12 ± 0.07 MeV and M B * - M B = 45.78 ± 0.35 
MeV. Both data may be used to obtain M H . c -M Scc and M s * b -M Sbb from Eq. J35J). If Q ^ Q', 
we use at NLL accuracy calculated in 3^|, and m b = M T (is)/2 and m c = Mj/^/2. For 
M=* — Mc we obtain about 107 MeV from the D data and about 133 MeV from the B 
data. Taking the average and estimating Aqcd/^c ~ oi1{m c a s ) ~ 0.3, our result is: 

M~* c - M Ecc = 120 ± 40 MeV. (36) 

Similarly for M^* b — M^ bb we obtain about 27 MeV from the D data and about 34 MeV from 
the B data. Taking only the estimate based on the B data, because affected by the smaller 
uncertainty Aqcd/^6 ~ «s( m b a s) ~ 0.1, our result is 

Af=. - Mn., = 34 ± 4 MeV. (37) 

^bb ^bb \ I 

These results compare well with the quenched QCD lattice simulation of jl^j], whose result 
is M=. - M= = 89 ± 15 MeV, and of % , whose result is M=. - Af= = 80 ± 10 + ? MeV, and 
with the quenched NRQCD lattice simulations of [8J] and (U|, whose results for bbq baryons 
are _M H * t -.M-,, = 20 ± 6+3 MeV and M B * - M S6ij = 20 ± 6+^ MeV respectively. The figures 



of 



<| and [lj 



refer to the lattice calculations at largest (3. 



IV. PNRQCD FOR QQQ BARYONS 

In this section we consider baryons formed by three heavy quarks of which at least two 
with the same mass mi = m 2 = m. Baryons of this type may be composed by bbb, bbc, bcc 



3 From private communications we know, however, that at least Tom Mehen and the authors of |2| were 
aware of it. 
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or ccc quarks. We define 



rriR = 2 m + m 3 , 



m 



2mm 3 



(38) 



m p = -, 



m x = 



R = 



m(xi + x 2 ) + m 3 x 3 



p = xi - x 2 , A = 



m R 
Xi + x 2 



(39) 



2 



- x 3 . 



m R 



There are, in principle, several physical scales that may play an important role in the dy- 
namics: the masses m R , m p and m x , which we assume to be of the same order, the typical 
relative three momenta of the heavy quarks, the typical kinetic energies and the scale of 
non-perturbative physics Aqcd- In the following, we will keep the discussion as simple as 
possible by not exploiting any possible hierarchy among the relative momenta and the kinetic 
energies. We will assume that the typical relative momenta of the heavy quarks, generically 
denoted by mv, are all much smaller than the heavy-quark masses and much larger than the 
kinetic energies, generically denoted by mv 2 . We may distinguish two situations. 

(A) The typical relative momenta of the heavy quarks are much larger than Aqcd- We 
call this situation weakly coupled. 

(B) The typical relative momenta of the heavy quarks are of the order of Aq CD . We call 
this situation strongly coupled. 

A. pNRQCD for weakly-coupled QQQ baryons 

1. Lagrangian and Degrees of Freedom 

If we assume that the typical distances p and A in the baryon, which are of order l/(mv), 
are much smaller than 1/Aqcd, then gluons of momentum or energy of order mv may be 
integrated out from NRQCD order by order in a s . The resulting EFT has light quarks, 
gluons of energy and momentum lower than mv (ultrasoft gluons), and heavy quarks as 
degrees of freedom. Gluons appearing in vertices involving heavy-quark fields are multipole 
expanded in p and A to ensure that they are ultrasoft. This corresponds to expanding in 
pAqcd ~ Aqcd/(wiO and AAqcd ~ Aqcd/(^^), or p{mv 2 ) ~ v and A (mv 2 ) ~ v. Like in 
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the case of doubly heavy baryons the number of allowed operators is consistently reduced if 
we choose to have a manifestly gauge-invariant Lagrangian. This may be done by projecting 
the Lagrangian on the heavy-heavy-heavy sector of the Fock space, by splitting the heavy- 
heavy- heavy fields into a singlet, two octet and a decuplet component, 

8 

Qu(x l ,t)Q 2j (x 2 , t)Q 3k (x 3 , t) ~ S(p, A, R, t) S ijk + Aa (p, A, R, t) o£j? 

a=l 

8 10 

+ Sa (p, A, R, t) Of/ fc + A V A, R, t) Aj fcl i,j = 1, 2, 3 , (40) 

a=l 6=1 

and by building the Lagrangian from these operators. The tensors S^, 0^£, Of and A_f - fe 

a 



are defined in appendix IB 21 S^- fe and A f - A are real and respectively totally antisymmetric 
and symmetric. We chose OJ^ and Of^ to be respectively antisymmetric and symmetric in 
the first two indices. The Lagrangian is constrained to satisfy all the symmetries of QCD. 
In particular, in the case m\ = m 2 that we consider here, it must be invariant under the 
exchange of the heavy quarks labeled 1 and 2. Under such transformation, p goes into 
— p and A goes into A. The gluon fields are even, because multipole expanded around the 
centre-of-mass of the heavy-heavy-heavy system. For what concerns the heavy-quark fields, 
from Eq. (|30|) it follows that S and Aa are even, because and 0^£ are odd under 
exchange i «-> j, and s a and A s are odd, because Of^ and Af jfc are even under exchange 
i <-> j. It is also useful to consider the combination of the above transformation with parity, 
(1 h 2) x P, which is also a symmetry of the Lagrangian. Under this transformation p goes 
into p and A goes into —A. The gluon fields transform like A^t, R) — > —R), which 
means that, up to reflection of the internal spatial coordinates, chromoelectric fields are odd 
and chromomagnetic fields are even. Up to reflection of the internal spatial coordinates, the 
heavy-quark fields transform like in the case of the (1 «-> 2) exchange. 

The resulting Lagrangian £ p nrqcd = £ p nrqcd(-R, t) at 0(X, p) in the multipole expansion 
(we also display at 0(l/m) the kinetic energy terms) is 

r -r 4. r 4- / ,(0,0) < 4- (au 

MpNRQCD — ^gluon T Might > ^pNRQCD "r ^pNRQCD "r ^pNRQCD ' \^ L ) 
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with £ g iuon and £n g ht defined in Eqs. (JHJ and (0) respectively and 



-(0,0) 
^pNRQCD 



+ o s t 



d 3 pd 3 X s f 



id - K 



(0)' 



S + At 



iD - V$ 



O' 



iD - V$ 



s + A f 



(0) 



A. 



V A p-EO s 



a=l v 
_ 8 _ / fabc I q Jafcc \ 

£ ( ' ; j [0 A ^p-gE b O Sc + O s ^p-gE b O Ac ] 

abc=l ' 
8 10 T / 3 \ 

EE E ^^fk )0^p.gE b A« 

ab=l 5=1 L \ii'jj'k=l / 

/ 3 \ 



-( E A^^^y^'fe) A*p-g&0 Aa 

\ii'jj'k=l I 



8 1 

- V^ eqA -7= [ Six ■ 9 Ea ° Aa + ° Aa] A • 9 Ea S] 

a=l * 

- v^i y (V 6c 2m ~ m3 + — 1 o a ^a • ^o A< 

abc=l ' 



aftc=l 
8 10 

+ t/ (0>1) V V 

afe=l ,5=1 



V fir &c 5m3 ~ 2m + — 1 Sa ^X-gE b O Sc 
f— ' V 6m R 2 / 



6m fi 

3 



77= E e ^ 

v ii'jj'k=l 



- V 



(0,1) 



~ ( T7f E ^!jk T ii' T ^i'j'k 

\ ii'jj'k=l / 

2m - 2m 3 £ £ / £ \ A 5t A • p£? a A 5 ' , 

0=1 <5<5'=l Wjjfc=l / 



A 5t A-^ b O Sa 



AA-EA mfi 



(1.0) 

pNRQCD 



d 3 pd 3 \ 5 1 " 



2m/j 2m p 2mA 



5 + At 



2mfl 2m p 2mA 



+ o s t 



+ 



2rriR 2m p 2mA 



s + A f 



+ 



2wifl 2m„ 2mA 



A 
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where the gauge fields in the covariant derivatives acting on the octets, A = (0 A1 , A2 , . . . , 
AS ) and s = (0 S1 , s2 , s8 ), and the decuplet, A = (A 1 , A 2 ,..., A 10 ), are under- 
stood in the octet and decuplet representation respectively. The dots in the last line of 
Eq. stand for terms that appear at orders higher than tree level and other 1/m terms, 
similar to those discussed for the doubly heavy baryon case. These terms are suppressed in 
the power counting with respect to the kinetic energy and the terms shown in Eqs. (j42j) and 

(JUS). 

The functions V are the Wilson coefficients of pNRQCD. They encode the contributions 
coming from gluons of energy or momentum of order mv . They are non- analytic functions 
of p and A. As we will discuss in the next section, at tree level we have 



V. 



(o,i) 



V, 



(0,1) 



Sp-EO s v O A pEO s 



V, 



(0,1) 



A p-EA 



sXeo a 

,(0) T/ (0) T/ (0) 



V 



(0,1) 



V 



(0,1) 



V 



(0,1) 



V 



(0,1) 



o a X-eo a o s Xeo s o s X-ea aX-ea 



1 



(45) 



while Vp, V^ U A ; , Vg a ' and get the first non- vanishing contribution at order a s . The coef- 
ficients in front of the operators D 2 R , V 2 , and V 2 are equal to 1, due to Poincare invariance 
or dynamical considerations similar to those developed in 31 1. 

The power counting of the Lagrangian (|4*T|) in the centre-of-mass frame goes as fol- 
lows: V A ,V P ~ mv, p, X ~ l/{mv), D R ~ A QC d,^ 2 , V s°o a < s ,a ~ mt>2 and E Br ^ 
Aq CD , {mv 2 ) 2 . The pNRQCD Lagrangian at leading order reads: 



^pNRQCD 



+o A t 

+ At 



d 3 pd 3 X \ S f 



v 2 V 2 X 

id + tt^- + 



2m r 



Vq 

2m x s 



(o) 



S 



iDr 



2m 



2m > 



V, 



(o) 



A + Sf 



X 



2m n 2m 



V, 



(o) 



A 



0> 



2m p 2m\ 



V, 



(0) 



f=i 



(46) 



0=1 



21 



2. Matching 



The matching from NRQCD to pNRQCD is performed by calculating Green functions 
in the two theories and imposing that they are equal order by order in the inverse of the 
mass and in the multipole expansion. Since we are working here in the situation where 
the typical momentum transfer between the heavy quarks is larger than Aqcd ; we can in 
addition perform the matching order by order in a s . The procedure is analogous to that one 
discussed previously for the doubly heavy baryon case, which we follow closely. 

The matching at tree level may be performed by projecting the NRQCD Hamiltonian on 
the three-quark Fock space spanned by 

3 

d 3 Xl d 3 x 2 d 3 x 3 ^ f lQ2Q .M,x 2 ,x 3 )Q\\ Xl )Qi\x 2 )Q^(x 3 )\0), (47) 

ijk=l 

where $QiQ 2 <33( £C 1) x i, x ~i) is a 3 ® 3 £§) 3 tensor in colour space and a 2 <g> 2 <g> 2 tensor in 
spin space. After projection, all gluon fields are multipole expanded in p and A. In order to 
make gauge invariance explicit at the Lagrangian level, we decompose the three quark fields 
into a field S(p, A, R, t), which transforms like a colour singlet, two fields A (p, A, R, t) and 
s (p, A, R, t), which transform like octets, and a field A(p, A, R, t), which transforms like a 
decuplet: 

3 

$oJQaQs( a! l> aj 2> a! 3>*) = < Pii'{ x ^ R ^ t ) ( t ) jj'{ x 2^R]t)(j) kkl {x^R]t) 

i'j'k'=l 

' 8 

S(p, A, R, t) S iYfc , + ° Aa (P> A > *) &?v 

^ a=l 

8 10 \ 

+ Sa (p, A, R, t) O^y + AS (P> A ' *) &ii>» '( 48 ) 

o=l 5=1 / 

where S ijfe , O^, Of^ and A^- fe have been defined in appendix IB 21 and the Wilson string 
in Eq. (|15j). After projecting on (|48p the Lagrangian (|42|) - (j44|) with the matching conditions 
(1431) follows. 

The perturbative matching of the static potentials Vg°\ VqI, VqJ and V^ '* goes as follows. 
In NRQCD we compute static Green functions, whose initial and final states overlap with 
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the singlet, octet and decuplet fields in pNRQCD. A possible choice, working in a non-gauge 
invariant framework, is 

3 

Im= E (0\M^ k Q l (R,x u T/2)Q J (R,x 2 ,T/2)Q k (R,x 3 ,T/2) 

ijki'j'k'=l 

xM:;, Ql(R, Vl , -T/2)Q],(R, y 2 , -T/2)Qj,(H, y 3 , -T/2)|0), (49) 



(1) if m = s, iz = h, Mi Jk = % k , 

(2) if M = A , I% = I&, M^ = 0^, u,v = l,2,...,8, 

(3) if M = O s , I% = I%1, M? jk = Ol u k , u,« = l,2,...,8, 

(4) if M = A, I% = I% V , M^ = A^. fc , u,v = l,2,...,10, 

where Q(R,x,t) has been defined in Eq. (J23j) . Integrating out the heavy-quark fields from 
we obtain 

IZ = 8\ Xl - Vl )5\x 2 - y 2 )5\x 3 - y 3 )(0\(W^ QQ r\0) , (50) 
with Wqqq diagrammatically represented in Fig. El and explicitly given by 

3 

(W^ QQ r = P J2 4>w(R, x h T/2)(j) Vrl {T/2, -T/2, xx)(j) r/r ( Xl , R, -T/2) 

ijki'j'k'rstr's't'=l 

(51) 

x ^(R, x 2 , T/2)0 iV (T/2, -T/2, x 2 )4> s , s (x 2 , R, -T/2) 

x <f> kk ,{R, x 3 , T/2)<t> k , t ,{T/2, -T/2, x 3 )<j) t , t (x 3 , R, -T/2) M^ t . 

In the large T limit, the Green functions Is, Iq\, Iq V s and I™ are reduced to the singlet, octet 
and decuplet propagators of pNRQCD respectively. If we neglect subleading loop corrections 
to the pNRQCD side of the matching, we obtain: 

lim (OKW&gHO) = lim Z M (p, X) exp (-ivff (p, X) t) 

1 — >oo 1 — »oo \ / 

3 

x (0£M^, <MT/2, -T/2, R)<f> jf (T/2, -T/2, R)<j> kV {T/2, -T/2, R) M^|0>, (52) 

ijki'j'k'=l 
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FIG. 2: Static Wilson loop with edges x\ = (x\,T/2), x 2 = (x 2 ,T/2), x 3 = (x 3 ,T/2), yi = 
(xx, — T/2), y 2 = (x 2 , — T/2), y 3 = (x 3 , — T/2) and insertions of the tensors A4^ fc and At^*, fc/ in 
X = (R, T/2) and Y = (R, -T/2) respectively. 



where Zm is a normalization factor. At order a s , the result is 

1 1 1 





(P,A) = 


2 

— a s 
3 




(P, ^ = 


2 

Q!s 

3 


v$ 


(P, ^ = 


«s / 

3 VI 




(P, ^ = 


«s / 

3 VI 



|p| |A + p/2| |A — p/2| 

111 11 



\p\ 8|A + p/2| 8|A-p/2| 
5 1 5 1 



(53) 
(54) 
(55) 

11 1 \ 

+ p^-TT^T + T\ T7^ • (56) 



|p| 4|A + p/2| 4|A-p/2 
1 1 1 

R + |A + p/2| + |A-p/2| 

B. pNRQCD for strongly-coupled QQQ baryons 

1 . Lagrangian and Degrees of Freedom 



In the situation in which the typical distances p and A in the baryon are of the order 
1/Aqcd, the matching from NRQCD to pNRQCD cannot rely on perturbation theory any- 
more. Also, it is more difficult to identify the effective degrees of freedom of pNRQCD. 
Despite these difficulties, the situation appears pretty much similar to that one described for 



strongly-coupled quarkonium in 
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From the available lattice simulations (e.g. [19|, see 
Fig. El), it appears that the gluonic excitations between three static quarks develop an energy 
gap of about 1 GeV > Aqcd with respect to the lowest static energy. This means that all 
gluonic excitations between the heavy quarks are integrated out once we go to pNRQCD. 
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L min [ fm ] 



1.5 



FIG. 3: Lattice measurements of the three-quark static energies of the lowest state, Eq*\ and of 
the first gluonic excitation, E^°\ as a function of L m i n , the minimal total length of the flux tubes 



linking the three quarks. From 



pNRQCD in its simplest formulation, i.e. without light quark degrees of freedom, is, there- 
fore, as simple as a potential model. 4 It has the three-quark singlet field S(p, A, R, t) as the 
only degree of freedom and is described by a Lagrangian £ p nrqcd = £ p nrqcd(-R, t), which 
reads: 



£pNRQCD — J d 3 pd 3 \ 



V 2 V V 2 

2mjj 2m n 2m\ 



S. (57) 



The potential Vs may be organized in an expansion (not necessarily analytic [37]) in the 
inverse of the heavy-quark masses. In the following, we will consider the matching of the 
1/m potential, which, in the non-perturbative regime, may, in principle, be of the same order 
as the static potential, and the 1/m 2 spin- dependent potentials. 



4 The relevance of the energy gap in relation to the success of the quark model has also been stressed in 
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2. Matching 



The non-perturbative matching goes as in the quarkonium case discussed in 2^| to 
which we refer for further details. Here we only list some results. 
The singlet static potential is given by 



V< 0) (p,A) = lim ^ln(0|Wg OQ |0>, 



T^ooT 



(58) 



where Wqqq is the singlet Wilson loop defined in Eq. (jHTj) and shown in Fig. El Lattice 
evaluations of may be found in 2^|. A plot is shown in Fig. El 
The order 1/m potential is given by 



, s t/ (1,1) T/ (1 ' 3) 
m m 3 



(59) 



with 



i " poo 

V^' 1 \p,X) = --Y, / dtt {(gE{x h t) ■ gE{x h 0))) S cA 
i=i ^° 

dtt((gE(x 3 ,t)-gE(x 3 ,0)))l QQQ 



i r°° 

2io 



(60) 
(61) 



where the double brackets stand for the gauge field average in the presence of a static Wilson 
loop of infinite time length: 



= lim L-L- wv^ / 62) 

«Oi(*i)O a (ta)»JgoQ = (Oi{h)0 2 {t 2 ))) s QQQ - ((Oi(ti)))§ QQ ((02(t2)})5 Q Q (63) 

, T T 

with — > t\ > to > . 

2 ~ 2 



As in the quarkonium case [2^J, in the non-perturbative regime the 1/m potential may, in 
principle, be of order mv 2 , and therefore as important as the static potential. There are no 
available lattice data for this quantity. 

For the potentials responsible for the spin splittings of the heavy baryons, we obtain at 
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order 1/m 2 : 



(2, spin dep.) 
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fe/=i 
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dt '£((gB k (x i ,t)gB l (x i ,,0))) 



-iV 



c,QQQ"k "l 



fe/=l 



E yQiQi' + d Q iQil iT^T^QQo) ^ • ^^(^ - <*0 , (64) 



i>i'=l 



where T a ^T a ^ s > stands for two colour matrices T a inserted at the same time in the Wilson 
lines of spatial coordinates a;, and av respectively, and the matching coefficients dg.q 7 and 
cJq.q., have been calculated in appendix El The above expressions give at order ct s the 
well-known one-gluon exchange results |22|. The spin-dependent potentials have not been 
calculated on the lattice yet, differently from the quarkonium case, where such calculations 
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42]. 



have instead a long history |18|]. Model dependent predictions may be found in 
It is expected that these potentials satisfy some exact relations due to Poincare invariance 



of the type studied in 



31 



43, 



44j for the quarkonium case. 



V. CONCLUSION AND OUTLOOK 

This work is a first step in the direction of a complete study of baryons made of two or 
three heavy quarks in the framework of non-relativistic EFTs of QCD. For both types of 
baryons, we identify the degrees of freedom and write the pNRQCD Lagrangian appropriate 
to describe the system in the heavy-quark sector. In the doubly heavy baryon case this 
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represents an update of Ref . [13j , which, however, used a HQET framework. In the case of 
baryons made of three heavy quarks, we also provide non-perturbative expressions for some 
of the potentials. Relevantly for both types of systems, we calculate the one-loop matching 
of the 4-quark operators of lowest dimensionality. 

Several further developments are possible. For doubly heavy baryons, where data are 
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already available, an important step forward would consist in providing pNRQCD with a 
light-quark sector that fully implements chiral symmetry and chiral symmetry breaking ef- 
fects. One could then study, for instance, isospin splittings and transitions and also address 
a variety of decay and production processes. The pursuit of such a program of phenomeno- 
logical studies will, however, very much depend on the future of the experimental searches 
for these states. 

For what concerns heavy baryons made of three heavy quarks, in the absence of a dis- 
covery, lattice studies will remain the main source of information. First of all, it will be 
important to have at least the one-loop expressions for the heavy-baryon static potentials 
Vg°\ VqI, Vqs and (also for Vj^ and V^). This may lead to a precise comparison of 
short-range lattice data with perturbative QCD in the heavy-baryon sector. At three loop, 
the heavy-baryon static potentials exhibit an ultrasoft running like in the heavy- quarkonium 

(0) 

case [45]. The ultrasoft running of the singlet static potential Vg comes from the coupling 
with the octets A and s . The leading logarithmic contribution at order a* is 

/V (0) - T/ (0)X2 

svP = t*\'(v<»-v$>A\> oA 



s " 9tT vo a r s j 47r/i2 



where /i is the ultrasoft factorization scale. 

Let us comment on the renormalon singularities affecting the perturbative series of the 
baryonic static potentials. These must cancel in physical observables. In the quarkonium 
case, the renormalon of the static potential cancels against twice the renormalon affecting 
the heavy-quark pole masses (see e.g. Q). From Eq. (j^ one can read that the order A QCD 
renormalon affecting is 3 x 1/2 that one of the static potential in the quarkonium case. 
Indeed, in the expression of the baryon mass it cancels against three times the renormalon 
affecting the heavy-quark pole masses. Similarly, in the doubly heavy baryon case, from 
Eq. (}2Tj) we have that the renormalon of order Aqcd affecting is 1/2 that one of the 
static potential in the quarkonium case. In the expression of the baryon mass, it cancels 
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against the renormalon affecting the A parameter of the HQET and the two heavy-quark 
masses. 

Concerning the energies of gluonic excitations from three static sources, in the short-range 
they are expected to behave like the singlet potential if they are singlet plus glueball 
states, or like the octet or decuplet potentials (|54 )1 -(|56 |1 if they are hybrid states. Only if E^ 
corresponds to the first case, the Coulomb contribution is expected to cancel in — E^\ 
which is the difference between the energy of the first excited state and the ground state. 
This could be in contradiction with a statement in Ref. 3], where the Coulomb contribution 
is said to cancel in the difference without any further specification. Like in the quarkonium 
case [29] , it is expected that the ordering of the levels of the gluonic excitations in the short 
range is dictated by the correlation lengths of some gluonic operators. If we assume that 
correlation lengths of operators of higher dimensions are suppressed and if we consider that 
there is a singlet channel only in 8® 8 but not in 10® 8, then, in the short range, the leading 
gluonic excitation is expected to come from the coupling of an octet heavy-quark state with 
a gluon field. It would be interesting to investigate if the first gluonic excitation shown in 
Fig. El is such an octet hybrid, and in this case what kind of octet. If it is not an octet 
hybrid, then likely it exists a lower gluonic excitation that still needs to be identified. 

Finally, in the perspective of a future spectroscopy of baryons made of three heavy quarks, 
it may become important to have a lattice determination of the spin-dependent potentials. 
Moreover, as the history of quarkonium suggests, spin-dependent observables may provide 
an excellent insight into the quark- confinement mechanism in the baryonic sector. 
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APPENDIX A: 1-LOOP MATCHING OF 4- QUARK OPERATORS OF DIMEN- 
SION 6 



The only graphs contributing to the 1-loop matching of the 4-quark operators of dimension 



6 are displayed in 
masses studied in 



ig. |3J The situation is similar to the quark-antiquark case with different 



47|. 





FIG. 4: Feynman diagrams contributing to the 1-loop matching of the 4-quark operators of dimen- 
sion 6. 



In the case of 2 different quarks of masses rrih and m^ (h ^ h!) we obtain in the MS 
scheme: 
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where C4 = N c — 3 and CV = (A*" 2 — 1)/(2N C ) = 4/3. For m/j = = m the above formulas 
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become 



dQ Q = M^-^)«s(m^-jf) , (A5) 



dQQ = C F [^-C F )ai, (A6) 



2 2\ 1 „ / m 2 23 



Cao; 2 (. m 2 



<*3b = 2C W - In— + 7 . (AS 



Working in D dimensions, we have used the prescription eij^ijk = (D — 1)(D — 2)(D — 3). 
If the prescription e^e^t = (D — 1)(D — 2) of j^J is instead used, this amounts to changing 

d Q h Q h > d Q h Q h > + C^/ 2 - 

APPENDIX B: GROUP FACTORS 
1. Multiplet Tensors: 3 ® 3 

The product of two triplet representations of SU (3) may be decomposed into the sum of 
an antitriplet and a sextet representation: 3 ® 3 = 3 © 6. A possible matrix representation 
for the antitriplet (T^, i, i,j = 1, 2, 3) and the sextet (E?-, a = 1, 2, 6 and i,j = 1, 2, 3) is 

Xii = -y| e «i' ( B1 ) 

Ull = = £33 = 1, 

S2 x~*2 — ^ /'T3 r ^ 

—12 — ^21 — ^13 — ^31 — ^23 — ^32 ~ 

all other entries are zero. 

Both Tjj and are real; T|- is totally antisymmetric and S^- totally symmetric. They 
satisfy the orthogonality and normalization relations: 

ET^Tj = 5«', ££^ = ^', ^ = 0. (B3) 

ij=l ij=l ij=l 
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2. Multiplet Tensors: 3 <g> 3 <g> 3 



The product of three triplet representations of SU (3) may be decomposed into the sum of 
a singlet, two octet and a decuplet representation: 3®3(g>3 = 1©8©8©10. A possible matrix 
representation for the singlet (S ijk , k = 1,2, 3), the octets (Q^% and Ofj^., a = 1,2, 8, 
i,j, k = 1,2, 3) and the decuplet (A*. fc , 5 = 1,2, 10, i,j, k = 1,2, 3) is 



1 3 

— o ^ e «jnAfc n , (B5) 
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1 3 
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all other entries are zero, 



(B7) 



where A a are the Gell-Mann matrices. S ijk and A.f jk are real; S ijk is totally antisymmetric 
and A.f jk totally symmetric. The octets and 0^ a fc have been chosen to be respectively 
antisymmetric and symmetric in the first two indices. The matrices satisfy the orthogonality 
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and normalization relations: 

3 3 3 3 

EC C -| \ riAo* 0^ a ' X aa ' \ r nSa* riSa' A aa ' \ r A <5 A 5' r<5<5' 

&ijk — ijk — L ' ±±-ijk ±Lijk — i 2-^i — l i k — i i k ~ ' 2.^1 — { j k — l 3 k ~ ' 

ijk=l ijk=l ijk=l ijk=l 

3 3 3 

SjjfcO^£ = ^ S ijA . = ^ Sijk^ijk = > (B8) 

ijk=l ijk=l ijk=l 

3 3 3 

E/-vAa*pjSa' \ r r\Aa * a <5 \ r z-^S a * a <5 r> 

k^jfc — —ijk ==ijk — 2.^1 —V k — l 3 k ~ ' 

ijk=l ijk=l ijk=l 

3. SU(3) Representations 

Here we list our choice of matrix representations for the SU(3) generators in the 3, 3, 6, 
8, 10 representations: 

w 

rpa — rpa — 

- 3 - 2 ' 

\aT 

3 2 

3 

(rpa\ — \ v" 7 \ a V" 7 ' 
l J 6V = 2^ =^ij A jk±±kii 

ijk=l 

(T 8 a ) bc = zf bac , 
3 3 

ii'jk=l 







(B9) 






(BIO) 


o,o' = \,2,. 


..,6, 


(Bll) 


6,c=l,2,.. 


.,8, 


(B12) 


5,5' = 1,2,.. 


.,10, 


(B13) 



where a = 1, 2, . . . , 8. 



APPENDIX C: SPIN STATES 



Let us consider a meson made by a heavy antiquark Q and a light quark g. We denote by 
\P*; S Z Q q ) the lowest Sg q = 1 states (SI = 1, 0, —1), and by |P; 0) the lowest = state. 
The heavy antiquark content of the states may be made explicit by writing: 

|P*;1> = J d 3 RQl(R)\St = l/2), (CI) 

|P*;0> = J d*R ^(Ql(R)\St = -1/2) + Ql(R)\Sl = 1/2)), (C2) 
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|P*; -1) = J d 3 R Ql(R)\St = -1/2), (C3) 

|P;0) = Jd 3 R-j=(Ql(R)\St = -1/2) - Ql{R)\Sf = 1/2)). (C4) 

In the case of the lowest doubly heavy baryon states, we denote by |S*; Sqq ? ) the Sqq q = 

3/2 states {S z QQq = ±3/2, ±1/2), and by |S;±l/2) the S QQq = 1/2 states. The heavy 
antitriplet content of the states may be made explicit by writing: 

|S*;3/2) = J d 3 Rd 3 r VQQ (r)Tl(r,R)\S* = 1/2), (C5) 
IS*; 1/2) = J d 3 Rd 3 r^ QQ {r) Uh[« (r, R) |Sf = -1/2) 

±y / |T t (r, J R)|Sf = l/2)Y (C6) 

|S*;-l/2) = J d 3 Rd 3 r p QQ (r) (^TS(r,R)\Sf = -l/2) 



where 



+ v /V(r, J R)|Sr = l/2)j, (C7) 

-3/2) = Jd 3 Rd 3 r VQQ (r)Tl(r,R)\St = -1/2), (Cf 
5; 1/2) = | d 3 Pd 3 r^ (r) ^4 ( r , H)|Sf = -1/2) 

|3j(r,JR)|S? = 1/2)Y (C9) 

|S;-l/2) = J d 3 Rd 3 r^ QQ {r) (J^T^r, R)\Sf = -1/2) 

-y|2l(r,fl)|Sf = l/2)V (CIO) 

| rf 3 r^ Q (rV QQ (r) = l. (Cll) 
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